Melting Pattern of Diquark Condensates in Quark Matter 
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Thermal color superconducting phase transitions in high density three-flavor quark matter are 
investigated in the Ginzburg-Landau approach. Effects of nonzero strange quark mass, electric and 
color charge neutrality, and direct instantons are considered. Weak coupling calculations show that 
an interplay between the mass and electric neutrality effects near the critical temperature gives rise 
to three successive second-order phase transitions as the temperature increases: a modified color- 
flavor locked (mCFL) phase (ud, ds, and us pairings) — > a "dSC" phase (ud and ds pairings) — > 
an isoscalar pairing phase (ud pairing) — > a normal phase (no pairing). The dSC phase is novel in 
the sense that while all eight gluons are Meissner screened as in the mCFL phase, three out of nine 
quark quasiparticles are gapless. 

PACS numbers: 12.38.-t,12.38.Mh,26.60.+c 



Unraveling the phase structure at high baryon den- 
sity is one of the most challenging problems in quantum 
chromodynamics (QCD). Among others, color supercon- 
ductivity in cold dense quark matter has been discussed 
from various viewpoints [l|,|2|. In relation to real systems 
such as newly born compact stars in stellar collapse, it is 
important to study the color superconductivity not only 
as a function of the quark chemical potential \x but also 
as a function of the temperature T. This is because the 
possible presence of color superconducting quark matter 
in a star affects the star's thermal evolution |3|. 

The purpose of this Letter is to investigate phase tran- 
sitions in color superconducting quark matter with three 
flavors (uds) and three colors (RGB) near the transi- 
tion temperatures. We consider a realistic situation in 
which nonzero strange quark mass m s , electric and color 
charge neutrality, and direct instantons take effect. As we 
shall see in weak coupling (m s , Aqcd <C //), the effects of 
nonzero m s and electric neutrality are important in that 
they induce multiple phase transitions that change the 
pattern of diquark pairing as T increases. In particular, 
we find a new phase, which we call "dSC," as an interface 
between a modified type of color-flavor locked (mCFL) 
phase and an isoscalar two-flavor (2SC) phase. 

Throughout this Letter, we adopt the Ginzburg- 
Landau (GL) approach near the transition temperatures, 
which was previously used to study the massless three- 
flavor case and is a more advantageous framework 
to weak coupling calculations than other mean-field ap- 
proaches 13 • In a realistic situation, the GL potential 
acquires the following corrections. First of all, nonzero 
m s affects the potential through the s quark propagator 
H in such a way as to lower the temperature at which a 
diquark condensate with s quarks dissolves. This is be- 
cause the pairing interaction due to one-gluon exchange 
is effectively diminished by m s if the pair contains the s 
quark. Secondly, when quark matter with nonzero m s is 
beta equilibrated and neutralized by electrons near the 
transition temperatures, the chemical potentials between 



d, s quarks and u quarks differ. Through this chemi- 
cal potential difference, the GL potential acquires an- 
other m s dependence, which is the essential origin of the 
dSC phase. Thirdly, the instanton contribution gives an 
m 8 dependence to the GL potential through the effective 
four-fermion interaction proportional to m s . We finally 
note that in weak coupling, color neutrality makes neg- 
ligible difference in the GL potential near the transition 
temperatures Q. 

We assume that diquark pairing takes place in the 
color-flavor antisymmetric channel with J p — + , which 
is predicted to be the most attractive channel in weak 
coupling 0- In this case, the pairing gap of a quark 
of color b and flavor j with that of color c and fla- 
vor k has the form cj) bc j k = e Q 6 C eyfc(d a ) 1 Q{. Here 
the 3x 3 matrix (d a )' transforms as a vector under 
G = SU(3) C x SU(3) L+R x U(1) B and belongs to the 
(3*, 3*) representation of SU(3) C x SU(3) L+R . 

For a homogeneous system of massless quarks 
{fnu.d,s = 0), the GL potential is invariant under 
SU(3)c x SU(3)l+r and valid near the critical temper- 
ature, T c , common to all states belongin g to the channel 
considered here. This potential reads [4|, |l£| 

S = |d a | 2 + |d a | 2 ) 2 + ft £ |d a * • d,| 2 , (I ) 

a a ab 

where (d a y = (d%, d d a , d s a ) and the inner product is taken 
for flavor indices. In the weak coupling regime, the coef- 
ficients are 
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N(p) = p, 



W(n)t = a t, (2) 



where N(fi) — /i 2 /27r 2 is the density of states at the Fermi 
surface, and t = (T — T c )/T c is the reduced temperature. 
With the parameters J2J), one finds a second order phase 
transition at T = T c from the CFL phase (d l a oc 8 l a ) to 
the normal phase (d l a = 0) in mean-field theory j^j. 

Let us now consider the effect of a nonzero m s in 
the quark propagator on the GL potential. We assume 



2 



m Ul d = for simplicity and consider the high density 
regime, m s <C /i. Near T c the leading effect of m s is 
to modify the quadratic term in the GL potential IjTJl. 
The corrections to the quartic terms are subleading and 
negligible. Since m s affects only us and ds pairings, the 
correction to the quadratic term has the form 



(K\ 2 + \d d a \ 2 ) = eJ2(\da\ 2 -K\ 2 ). 



(3) 



Note that this correction induces an asymmetry in the 
flavor structure of the CFL phase. 

In weak coupling, e can be calculated by including m s 
in the Nambu-Gor'kov quark propagator when evaluat- 
ing the thermodynamic potential. Following Ref. p|, we 
expand the thermodynamic potential not only in d l a but 
also in m s up to 0(m 2 ), and obtain 
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Here the dimensionless parameter a is given by 
f3ir 2 
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(4) 
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where g is the QCD coupling constant. As long as a <C 1, 
which is relevant at asymptotically high density, the fol- 
lowing GL analysis near T c is valid. In the latter es- 
timate in Eq. (@J we use the leading-order result in g, 
ln(T c /^) 37r 2 /(\/2ff) B- This behavior of T c orig- 
inates from the long-range color magnetic interaction 
which prevails in the relativistic regime. As a result of 
the modification by m s to ln(T c //x), e has a positive value 
such that ud pairing is favored over us and ds pairings. 
Consequently, the CFL phase becomes asymmetric in fla- 
vor space and its critical temperature is lowered, leading 
to the appearance of the 2SC phase (d l a oc 5 ls ) just below 
T c [1]. Note that Eq. © has no effect on the 2SC phase. 
We also note that T c itself is modified by m s through the 
modification of the normal medium as T c (l + 0{ga)). 

We turn to discuss effects of charge neutrality which 
also depend on m s as mentioned above. Under beta equi- 
librium and charge neutrality, the electron chemical po- 
tential fi e and the shift Sfii of the chemical potential of 
flavor i from the average (/z) are related by <5/i,; = —qifi e , 
with electric charge We see from Ref. that the 
thermodynamic-potential correction due to 8/ii has the 
form 



^Ei d «i 2 -EKi 2 )- 



(6) 



In weak coupling, where one may regard normal quark 
matter and electrons as nonintcracting Fermi gases, jU e 
is related to m s as [i e = m 2 /A^L. This estimate is valid 
in the vicinity of T c where corrections to /x e by a finite 



pairing gap affect only the quartic terms in the GL poten- 
tial. By combining this estimate with the weak coupling 
expression for r\j [i e given in Ref. 0, we obtain 



Tr 



a a. 



(7) 



Since 77 > 0, ds pairing is more favorable than ud and us 
pairings. This feature stems from the modification by dfii 
to the exponential factor of T c , exp[— 3n 2 /(y/2g)], which 
tends to increase (decrease) the critical temperature for 
ij pairing when + 5[ij > 0(< 0). 

We consider color neutrality of the system as well. In 
contrast to the case at T = 0, however, it affects only 
the quartic terms in the GL potential through possi- 
ble chemical potential differences between colors 0, ^] , 
and in weak coupling its magnitude is suppressed by 
C((T c /g/i) 2 ) compared to the leading quartic terms. 
Thus color neutrality has no essential consequence to the 
phase transitions considered in this Letter. A major dif- 
ference between corrections of the charge neutrality and 
the color neutrality is that the former shifts the quark 
chemical potentials even in the normal phase, while the 
latter works only when the pairing occurs. This is why 
the former is more important than the latter near T c . 

The direct instanton at nonzero m s , which induces 
an effective four-fermion interaction between u and d 
quarks 0, leads to a quadratic term in the GL poten- 
tial, C J2 a Mai 2 - A n explicit weak coupling calculation 
shows that C ~ -ao(m s /^)(A QCD /^) 9 (l/5) 14 . The neg- 
ative sign indicates that the instanton favors ud pairing 
as does one-gluon exchange [see Eq. ©], but the magni- 
tude of C is highly suppressed at high densities. Hereafter 
we will thus ignore instanton effects. 

Since the two effects of nonzero m s , characterized by 
Eqs. © and ©, favor ud pairing and ds pairing, respec- 
tively, the finite temperature transition from the CFL 
to the normal phase at m s = is expected to be signifi- 
cantly modified. In fact, successive color-flavor unlocking 
transitions take place instead of a simultaneous unlocking 
of all color-flavor combinations. To describe this hierar- 
chical thermal unlocking, it is convenient to introduce a 
parameterization, 



(8) 



where Ai : 2,3 are assumed without loss of generality to be 
real. We also name the phases for later convenience as 




Ai,2,3 + 

Ai = 0, A 2 , 3 ± 

A 2 = 0, Ai, 3 + 

A!, 2 = 0, A 3 ^0 



mCFL, 
uSC, 
dSC, 
2SC, 



(9) 



where dSC (uSC) stands for superconductivity in which 
for d (u) quarks all three colors are involved in the pair- 
ing. 
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FIG. 1: Transition temperatures of the three-flavor color su- 
perconductor in weak coupling: (a) all quarks are massless; 
(b) nonzero m s in the quark propagator is considered; (c) 
electric charge neutrality is further imposed. The numbers 
attached to the arrows are in units of <jT c . 



In terms of the parameterization Q8J, the GL potential 
with corrections of 0(m 2 s ) to the quadratic term, Eqs. 
and iJBJ, reads 



+ f3 1 (A 2 1 + Al + Al) 2 
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where a' = a + e + 3 . 

We proceed to analyse the phase structure dictated 
by Eq. i fTUj l with the weak coupling parameters ©, Q, 
and J7J up to leading order in g. In Figs. 1 and 2 the 
results obtained by solving the coupled algebraic equa- 
tions, dS/dAi£,3 = 0, are summarized. Figure 1(a) 
shows the second-order phase transition, CFL — > nor- 
mal for m s = 0. Figures l(b,c) represent how the phase 
transitions and their critical temperatures bifurcate as 
we introduce (b) effects of a nonzero m s in the quark 
propagator and then (c) effects of charge neutrality. In 
case (b), two second-order phase transitions arise, mCFL 
(with Ai = A 2 ) -► 2SC at T = T c s = (1 - 4ct)T c , and 
2SC — > normal at T = T c . In case (c), there arises three 
successive second-order phase transitions, mCFL — > dSC 
at T = T l c , dSC -► 2SC at T = T c n , and 2SC -> normal 
at T = T c m . Shown in Fig. 2 is the T-dependence of 
the gaps Ai, 2 ,3 for case (c). All the gaps are continuous 
functions of T, but their slopes are discontinuous at the 
critical points, which reflects the second order nature of 
the transitions in the mean-field treatment of Eq. I|10(l . 

We may understand the bifurcation of the transition 
temperatures in weak coupling as follows. In the mass- 
less case (a), T c is degenerate between the CFL and 2SC 
phases, the chemical potential is common to all three fla- 
vors and colors, and the CFL phase is more favorable 
than the 2SC phase below T c . As one goes from (a) 
to (b), nonzero m s sets in, which tends to suppress the 
pairing interactions including the s quark. The critical 
temperature for the CFL phase is then lowered, which 
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FIG. 2: A schematic illustration of the gaps squared as a 
function of T. 



allows the 2SC phase to appear at temperatures between 
T c s and T c . As one goes from (b) to (c), charge neutrality 
sets in, which acts to decrease the chemical potential of u 
quarks by (2/3)/i e and increase that of d and s quarks by 
(l/3)/i e . Since fi e > 0, the average chemical potential of 
quarks involved in ds pairing increases, while those in ud 
and us pairings decrease equally. Accordingly, the tran- 
sition temperatures further change from T c to T™ and 
from T c s to T\ and 7; 11 . 

Now we examine in more detail how the color-flavor 
unlocking in case (c) proceeds with increasing T from 
the region below T l c . 

(i) Just below XjJ, we have a CFL-like phase, but the three 
gaps take different values, with an order A3 > Ai > 
A2 7^ (the mCFL phase) . The reason why this order is 
realized can be understood from the GL potential 11011 . 
The e-term and r^-term in Eq. (|ll)f) tend to destabilize us 
pairing (A2) relative to ud pairing (A3) and ds pairing 
(Ai). Since e > r](> 0), ds pairing is destabilized more 
effectively than ud pairing. The value of each gap in the 
mCFL phase reads 
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(11) 

(12) 
(13) 



The mCFL phase has only a global symmetry 
U(1)c+l+r x U(1)c+l+r in contrast to the global sym- 
metry SU(3)c+l+r in the CFL phase with m u ^d, s = 0. 
There are no gapless quark excitations in both mCFL and 
CFL phases. As T increases, the first unlocking transi- 
tion, the unlocking of A2 (the pairing between Bu and 
Rs quarks), takes place at the critical temperature, 



V i 1 - yo- 1 r;. 



(14) 
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(ii) For T l c <T < T c n , A 2 = and 
"o [T c -T 



A 2 

^3 



6/3 V T c 3 
QtQ (T c — T 7 
6/3 



(15) 
(16) 



In this phase, we have only ud and ds pairings (the dSC 
phase), and there is a manifest symmetry, U(\)c+l+r x 
U(1) c +l+r x U (1) c +v+b x U (l) c +v+B- By diagonaliz- 
ing the inverse quark propagator in the Nambu-Gor'kov 
representation, we find three gapless quark excitations in 
the color-flavor combinations: Bu, Rs, and a linear com- 
bination of Ru and Bs. At T = Tj? 1 , the second unlocking 
transition, the unlocking of Ai (the pairing between Gs 
and Bd quarks), takes place at the critical temperature, 



t" i i - ^ | r, 



(17) 



(iii) For T c n < T < T c m , one finds the 2SC phase, which 
has only ud pairing with 



2 _ ao ( T c -T _ 1 
~T, 3 CT 



A^ = — 

3 4/3 



(18) 



The 2SC phase has a symmetry SU(2) C x SU(2) L+R x 
t/(l)c+B x U(1)l+r+b @ In this phase the s quark 
and _B quark excitations are gapless. The final unlocking 
transition where A3 (the pairing between Rd and Gu 
quarks) vanishes occurs at 



T. 



in 



1 



1 



T c 



(19) 



Above TjJ 11 , the system is in the normal phase. 

In Table I, we summarize the symmetry and the gapless 
quark modes in each phase discussed above. The num- 
ber of gluons having nonzero Meissner masses, which is 
related to the remaining color symmetry, is also shown 
[lij . We note that more gapless quark modes may ap- 
pear if the system is in the close vicinity of T\, T} 1 , and 
T c m where A 2 , Ai, and A 3 are less than ~ m 2 Ju (l3|. 



TABLE I: The symmetry, the gapless quark modes, and the 
number of Meissner screened gluons in the mCFL, dSC, and 
2SC phases. The gapless quark mode (Ru, Bs) denotes the 
linear combination of Ru and Bs quarks. 





Symmetry 


Gapless 
quark modes 


Number of 
massive gluons 


mCFL 


[U(l)f 


none 


8 


dSC 


[11(1)]* 


Bu, Rs 
(Ru,Bs) 


8 


2SC 


[SU(2)f x [U(l)f 


Bu, Bd, Bs, 
Rs, Gs 


5 



So far, we have studied the phase transitions in the 
mean-field level. In weak coupling, as shown in Ref. [(| 
in the massless limit, thermally fluctuating gauge fields 
could change the order of the transitions described in 
Figs. 1 and 2. A detailed account on this effect will be 
reported elsewhere 14] ; here we recapitulate the impor- 
tant results. First, the second order transition, mCFL 
— > dSC, remains second order even in the presence of the 
thermal gluon fluctuations. This is because all eight glu- 
ons are Meissner screened at T = T\ and thus cannot 
induce a cubic term with respect to the order parameter 
in the GL potential. On the other hand, the transitions, 
dSC — > 2SC and 2SC — > normal, become weak first or- 
der since some gluons, which are massless in the high 
temperature phase, become Meissner screened in the low 
temperature phase (Table I). 

In summary, we have investigated color-flavor un- 
lockings at finite temperatures taking into account the 
strange quark mass and charge neutrality in the GL ap- 
proach. We find three successive unlocking transitions, 
mCFL — > dSC — > 2SC-^ normal, occurring in weak cou- 
pling. Most remarkably, the dSC phase appears between 
the mCFL and 2SC phases. In this phase all eight gluons 
are Meissner screened and the three quark excitations are 
gapless. The question of how the phase structure of neu- 
tral quark matter obtained near T c is connected to that 
at T = 0, which approaches the CFL phase in the limit 
of high density, is an interesting open problem. 
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